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Abstract. Electronic realizations of neurons are of great interest as building 
blocks for neuromorphic computation. Electronic neurons should send signals 
into the input and output lines when subject to an input signal exceeding a given 
threshold, in such a way that they may affect all other parts of a neural network. 
Here, we propose a design for a neuron that is based on molecular-electronics 
components and thus promises a very high level of integration. We employ the 
Monte Carlo technique to simulate typical time evolutions of this system and 
thereby show that it indeed functions as a neuron. 



PACS numbers: 73.63.-b, 84.35.+i, 85.35.-p, 73.23.Hk 
1. Introduction 

Neurons are the fundamental building blocks of information processing by the nervous 
system of every multicellular animal on Earth. They are cells of various degrees of 
sophistication serving the ultimate goal of processing and transferring information to 
adjacent neurons, and thus coordinating all of the animal's activities. 

Reproducing some, if not all, of the capabilities of a biological neural network 
using electronic components has always been a goal of research into artificial 
intelligence pQ. The benefits would range from massively parallel unconventional 
computing [2] to electronics that adapts to the type of signal that is fed into the 
network [3j 0], to name just a few. Traditional solid-state transistors can indeed 
mimic the main characteristics of a neuron (see, e.g., [5]). However, they suffer 
from integration limits and therefore are unlikely to reach the scalability of biological 
brains. In this respect, neurons based on molecular electronics appear to be a viable 
alternative. For instance, recent work has demonstrated an organic-nanoparticlc 
transistor that operates as a spiking neuron [5]. However, the design [B] relies on a 
combination of a thin film of pentacene molecules and gold nanoparticles. To simplify 
its fabrication and increase its integration capabilities it would be desirable to instead 
design an all-molecular circuit that operates as a neuron. Such a neuron could serve 
as a building block for highly- integrated neural computers [7J. 

In the present work we propose a type of molecular neuron based on the 
mechanism of Franck-Condon blockade [H El HH1 [HJ [12l [13l HU [HI HH] and show 
that it reproduces all the main features of a spiking neuron. Our guiding principles 
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are the following: (i) the design should be as simple as possible, and (ii) it should 
not require excessive fine tuning of parameters in order to function. To start with, we 
have to specify how we want the device to behave. Typical molecular-electronics 
components suggested in the literature pump electrons between source and drain 
electrodes, controlled by a gate voltage. Hence, the input signal for the artificial 
neuron will likely be a current. Furthermore, we require the neuron to fire sharp 
voltage spikes into its output line but also into its input line when the input current 
exceeds a certain threshold, so that it may affect all other parts of the network [3]. 
Below the threshold, the neuron should be quiescent. 

We thus need a molecular device that generates sharp spikes. Our central 
idea is to use a molecular transistor in which the electrons are strongly coupled 
to a vibrational mode. This coupling can lead to Franck-Condon (FC) blockade 
[HI El [13 dl H21 H31 H31 HH US]- The essential physics of this phenomenon is the 
following: The relevant vibrational mode is described by a normal coordinate x. 
The potential (deformation) energy Vd e f(aO is a different function of x for different 
charge states. In particular, for strong electron-vibron coupling, the value x = x m in 
minimizing Vdcf (x) shows a large shift between charge states. The molecule is in some 
vibrational state at time t, described by a wave function tp{x,t). Let us assume that 
it is in the vibrational ground state. Now, when an electron tunnels in or out, the 
potential energy Vdcf (x) is suddenly switched to a quite different function, with shifted 
minimum. Hence, the molecule finds itself in a state far from the vibrational ground 
state for the new potential. More precisely, the overlap between the old and new 
ground states is very small. The overlap between other low-energy eigenstates with 
respect to the two potentials is also suppressed. These overlaps are called FC matrix 
elements [T71 [HI QZ3 HH HI] • Due to the small FC matrix elements, electronic tunneling 
transitions that are energetically possible can nevertheless be strongly suppressed, in 
particular for low bias voltages. Consequently, the current through a molecular device 
with strong electron-vibron coupling is suppressed at low bias. This effect is called 
Franck-Condon blockade. 

The dynamics in the FC regime is also unusual: The electrons tunnel in avalanches 
separated by quiescent intervals [HI [S] ■ The reason is that while the probability for an 
electron to tunnel into the molecular orbital is very small, when it finally does so and 
then tunnels out again, the molecule likely returns to a higher vibrational level since 
the FC matrix element for such a process is larger than for a return into the ground 
state. But this makes it much more likely for another electron to tunnel. Thus we 
get an avalance, until the system finally drops back into the vibrational ground state. 
Our idea is to use these avalanches to generate spikes when the neuron is active. 

The remainder of this paper is organized as follows: In section [2j we introduce 
the molecular circuit and qualitatively discuss its behavior. The model and simulation 
technique are discussed in section [3] In section [4] we present results for the time- 
dependent behavior of the neuron. We summarize the paper in section [5] 

2. Molecular circuit 

The circuit diagram is shown in figure [T] The active part of the circuit is a single- 
molecule transistor in the FC regime (dotted rectangle labeled "FC" in figure [T]). 
Disregarding the transistor labeled "aux 2" for the time being, the input current ii n 
leads to a voltage drop across the resistor i?j n . This voltage is applied to the gate of the 
FC molecular transistor. Thereby, the input current is able to switch this transistor 
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Figure 1. Molecular neuron circuit based on FC blockade. The three dotted 
rectangles enclose single-molecule devices, as explained in the text. The black 
circle marked "dot" denotes a quantum dot separated from ground by a tunneling 
barrier. 



between OFF and ON states. The OFF state does not have any molecular transition 
energies between the electrochemical potentials of the source and drain electrodes so 
that sequential tunneling is thermally suppressed. On the other hand, in the ON state 
such transitions exist, but are suppressed by small FC matrix elements. In the ON 
state, there is thus a current from the source electrode kept at fixed voltage Vtuas 1 to 
the drain electrode. As noted above, the current flows in avalanche-like bursts [11 [9]. 

Any avalanche dumps a relatively large number of carriers within a short time 
into a quantum dot attached to the drain (circle labeled "dot" in figure [TJ . This 
quantum dot may be realized by a small metallic cluster or a large organic molecule. 
The carriers leave the dot via another tunneling barrier (labeled "barrier" in figure [T]). 
The charge on the dot is coupled to the gate of two auxiliary molecular transistors, 
labeled "aux 1" and "aux 2" in figure [l] They are not in the FC regime and will 
in the following be modeled by a single orbital with strong Coulomb repulsion but 
not coupled to vibrations. The first auxiliary device controls current flowing from a 
source electrode at fixed voltage Vbias 2 to a drain electrode connected to the output 
line. When the dot is charged due to the FC molecule being active, a current can flow. 
By choosing large tunneling rates to the source and drain contacts, the device can, in 
principle, act as an amplifier. The current flows through a resistor R ou t to the ground, 
creating a voltage drop, which we use as the output signal (in the simulation, we will 
assume that the output line is current free). During an avalanche, the auxiliary device 
is active, a current flows, and there is a nonzero voltage drop. The result is a voltage 
appearing only in short time intervals, as required. We introduce a capacitor C ou t to 
ground to dampen the sharp voltage spikes, which would otherwise be (5-function-like. 
In a real setup, intrinsic capacitances would always lead to such a broadening. 

The second auxiliary device ( "aux 2" in figure [T]) essentially works like the first 
one, except that it inserts a current into the input line, when the FC molecule is active. 
The current leads to a voltage drop across the input resistor R- ln . The voltage spikes 
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are dampened by the gate capacitor of the FC molecule. We find that an additional 
capacitor is not required. The connection of both the gate of the FC molecule and 
the drain of the second auxiliary molecule to the same input line of course leads to 
feedback. If the feedback is too strong, a current avalanche through the FC molecule 
could be choked off immediately. However, we will show below that this undesirable 
behavior can be avoided by a suitable choice of circuit parameters. 

3. Simulation method 

In the field of transport through molecular devices, the choice of the theoretical 
approach depends on details of the system: If the hybridization between molecular 
orbitals and states in the electrodes is relatively weak so that it can be treated 
perturbatively, master-equation approaches [HI [HI [13 HI EH1 [HI H21 [201 H3] are 
appropriate. The strengths of interactions, such as the Coulomb repulsion between 
electrons in the molecule and the electron-vibron coupling, may be large in these 
approaches. On the other hand, if interactions are weak, methods based on non- 
equilibrium Green functions can be employed (2TJ [22]. These approaches are able 
to treat the hybridizations exactly and are thus not limited to small hybridizations. 
For our purposes, however, Green-function methods are not suitable since our neuron 
circuit relies on the strong electron-vibron interaction in the FC molecule. 

3.1. Model 

In this work, we assume the hybridizations to be weak and treat them in leading-order 
perturbation theory, i.e., in the sequential-tunneling approximation. We also assume 
that dephasing is rapid so that off-diagonal components in the reduced density matrices 
(coherences) of the molecular transistors can be ignored [5T] . The sequential-tunneling 
rates are then well known [17 1 151 ITU 1 IT3 1 ITS] . 

We take all molecular transistors to contain a single relevant molecular orbital. 
The Hamiltonian of the FC molecule reads 



where 0% creates an electron of spin a in the molecular orbital and w is the creation 
operator of a harmonic vibrational mode. The energy epc of the molecular orbital 
is shifted by the electric potential Vpc- The potential is coupled to the voltages 
Vbiasi, Vdotj and Vin applied to the source, drain, and gate electrodes, respectively, 
through the capacitances of these contacts, as further discussed below. The voltages 
and capacitances are indicated in figure [2]ja). huj v in ([I]) is the energy quantum of 
the harmonic oscillator and the dimensionless constant A describes the strength of the 
electron-vibron coupling. 

The eigenstates of the isolated molecule are written as \nq), where the electronic 
state is denoted by n = 0, t, I, tl f° r the empty, spin-up, spin-down, and doubly 
occupied state, respectively. The vibrational states are enumerated by the harmonic- 
oscillator quantum number q — 0,1,2,... The eigenenergies read |16j 
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Figure 2. Details of the circuit diagram of the molecular neuron, indicating 
quantities used in the text, (a) FC molecule, dot, and tunneling barrier, (b) 
Auxiliary molecule 1. (c) Auxiliary molecule 2. 
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where = 0, 1, 2 is the number of electrons in the state \nq). 

The sequential-tunneling rates from state \n'q') to state \nq) involving an electron 
tunneling out of the molecule into the source or the drain electrode are 
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respectively, where T-pc is a bare tunneling rate determined by the hybridization and 
the density of states in the electrodes, f(x) is the Fermi function, D° n , = {n\d a \n') 
are matrix elements of the electronic annihilation operator, and 

\ Fqq ,\ 2 = ^ \^>-i<) e - x2 [L«>-s<(A 2 )] 2 (5) 

q>- 

are FC matrix elements squared [T71 [HI [TUJ [TH [TB]. Here, g< = min(q, q'), q > = 
max(g, q'), and L l Jx) are generalized Laguerre polynomials. In the interest of a simple 
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model, we assume the bare tunneling rate Tpc to be the same for source and drain. 
This assumption is not crucial for a functioning neuron. The corresponding rates for 
an electron tunneling into the molecule read 

Kt'^m = r FC f(E m - E n , q , + eV hias i) £ \D° n , n \ 2 \F q , q \ 2 (6) 

cr 

and 

R l X^ nq = Tfc f(E nq - E nW + eV dot ) J2 IK'nl 2 M. (7) 

(J 

For the tunneling barrier between the dot and ground we take the rates 

D Tbarricr 

^ cxp(eV/ dot /fc s T) - 1 W 
for an electron tunneling from the dot to ground, and 

D "Tb airier 

^ 1 - exp(-el/ dot /fc s T) W 

for the reverse process. Here, 7b a mcr is a bare tunneling rate. Note that this 
ansatz satisfies detailed balance and leads to ohmic behavior since the current is 

-^barrier — e_R_j. -)- eR^- — e '"/barrier I'dot ■ 

The Hamiltonians for the two auxiliary molecules describe a single orbital with 
strong Coulomb repulsion but no coupling to a vibrational mode, 

-ffaux u = ^{taux v - eV aux v ) d\d a + U aux v dfad^d^, (10) 

v = 1,2, with eigenenergies E* ux v = (e aux v - eV a ux v) n d + (<7 aux „/2) n d (n d - 1). The 
sequential-tunneling rates are of the same form as in ([3| , Q , ^ , and Q , with the 
FC matrix elements replaced by unity and the relevant parameters shown in figures 
|§b) andfgc). 

Finally, the on-site potentials at the positions of the three molecules and of the 
dot are obtained by solving the equations 

T/ _ CpCsHias 1 + CFC,d^dot + Cpc.g^m 

V FC - n r-ij n ' ) UJ-J 

^FC,s + L-FC,d + (-TCg 

T , _ C aux l, s Vbias 2 + C aux \^dVo\i\, + C aux l.gVdot /-. „n 

"aux 1 — ~ , — i — i j 

^aux 1,3 ~r ^aux l,d t ^aux l,g 

t t _ C a ux 2,s Vbias 2 + C aux 2,dMn + C aux 2,gVdot ( T 

Vaux2 — — , — — — — , (.loj 

L'aux 2,s T L/ aux 2,d T Wiix 2,g 

y _ Qdot + CFC.d^C + C a ux l,gV aux 1 + C aux 2, s V aux 2 

CFC.d + C aux i.g + C aux 2,g + Cbarrier 

where Qdot = — eiVdot is the charge on the dot. Note that we do not include the 
charges on the molecules explicitly in these equations because we have already done 
so in the Hamiltonians -ffpc i ^aux 1 j and -ff aux 2 • 



3.2. Monte Carlo simulations 

We are interested in the time-dependent behavior of the circuit. While the master 
equation can be used to study dynamics [HI HI E3 1231 [Ml ESI HH], it describes 
the dynamics of ensembles. We can gain more insight by considering typical time 
evolutions of a single system. For a molecular transistor in the FC regime, this has 
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been done in [HI El [10] . We obtain the time evolution by performing Monte Carlo 
simulations for the full circuit, using the transition rates given above. 

We employ real-time Monte Carlo simulations, which do not involve discretization 
of time. The state of the circuit is characterized by the state \nq) of the FC molecule, 
the states \n)\, \n)2 of the two auxiliary molecular transistors, and the occupation 
number -/Vdot of the dot. The harmonic-oscillator ladder is truncated at q max = 30, 
which does not affect the results since this value is never reached. At every step, we 
first calculate the transition rates from this state to all other possible states. The 
transition that actually happens is then selected pseudo-randomly with a probability 
given by its branching fraction, i.e., its rate divided by the sum of all rates (the total 
transition rate). The waiting time is drawn pseudo-randomly from an exponential 
distribution with mean given by the inverse of the total transition rate. 

To avoid inessential complications, we assume a constant incoming current I ln , 
which should be valid as long as this current is large compared to the additional 
current through the second auxiliary molecule. We also assume that the output line is 
current free, / ut = 0. Despite these simplifications, this is still a complicated model 
due to the time dependence introduced by the RC elements. For the output line we 
have Kirchhoff 's current law 

- eiV out = + Content, (15) 

-K ou t 

where iV ou t (t) is the total number of electrons inserted from the first auxiliary molecule 
into its drain electrode since an arbitrary reference time. N out (t) decreases if electrons 
flow into the device. If the time of the tunneling process is small compared to the 
other time scales of the circuit, the time derivative N ou t can be treated as a series of 
S functions. The solution of ( 15 ) is 

V out (t) = V m t(fl) e-*/^ - ~ e-*/ T °- / ' dt' e*'^°- N on t(t% (16) 

t^out Jo 

where r out = i?outC ou t is the characteristic time of the RC element. For the input 
line, the current law reads 

/in - eA> in = ^ + C FC ,g (Vin ~ V FC ) . (17) 
-Kin 

Here, the onsitc potential of the FC molecule, V F q, depends on time not only through 
the gate potential V- ln but also through the potential on the dot, Vdot- The effect of the 
fluctuating dot potential on V\ n via the FC molecule is smaller than the contribution 
from the current —eN- m and we neglect it for simplicity. We can then write 

fco«^|>. (is) 

where the partial derivative can be obtained in terms of the capacitances from the 
solution of ([TT])-([l4|. We thus find 

An - eA> in « ^ + C FC , g (l - ^) V in = ^ + C in i^, (19) 

-Kin V "'in / -Kin 

where we have defined an effective input capacitance C; n . The solution is then 

V in (t) = R m I m + [VLx(O) - i? in /in] e-*/ r » - — e-*/ r - / dt' e*'/ Ti " N in (l/), (20) 

^in Jo 

with r in = R\ n C;„. 
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The time dependences of Vi n and V on t imply that the transition rates change 
continuously in time and not only discretely at the time of tunneling transitions. It 
would be difficult to implement these time-dependent rates in the continuous-time 
approach. However, this is not necessary since the characteristic times T- ln and T out 
are long compared to the waiting time between tunneling events. Electrons are rapidly 
tunneling back and forth between the quantum dot and ground, with a characteristic 
time of l/7barricr- Since Ti n , r out 3> l/7barrier for our parameters, we can ignore the 
change of the rates during the Monte Carlo steps. 

Consequently, for a Monte Carlo step of duration At, the output voltage is 
updated according to 

V out (t + At) = V out {t) e ~ A */-o Ut _ e ANmU (21) 

L'out 

where AAT out G { — 1,0,1} is the number of electrons inserted into the drain in this 
step. Similarly, the input voltage is updated according to 

V in (t + At) = R in I m + [V in (t) - R in I in ] e- At ^ - ^- AN- m . (22) 



All the other voltages are then recalculated using (11)-(14| and the algorithm loops 
back to the calculations of the transition rates out of the new state. The Monte Carlo 
step is repeated until the simulated time t exceeds an equilibration time t; n ;t in order 
to get rid of any dependence on the initial state. Then, the simulated time is reset to 
zero and the state of the system and the various voltages are recorded at each Monte 
Carlo step until the simulated time exceeds t s im- 



4. Results and discussion 



A functional device is obtained for the following parameters: The energies in the 
molecular Hamiltonians are expressed in terms of a basic energy unit e as 

e FC = 12.9e , U FC = 500e , hu v = e , A = 4, (23) 
faux i = -3.6 e , £aux 2 = -5.5 e , f7 aux i = U aux 2 = 500 e , (24) 

where the values for the Coulomb repulsions prevent double occupancy and are thus 
effectively infinite. The temperature is taken to be ksT — 0.05 eo. The parameters for 
the FC molecule are thus essentially the same as used for the Monte Carlo simulations 
in [51 [3] ■ In particular, the electron- vibron coupling A = 4 is strong and the molecule 
should show FC blockade. The other parameters of the circuit are 
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= - Vbiasa = -3co/e, 


(25) 
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= o.oi r . 


(34) 
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Note that the contacts of each molecule to its source and drain electrodes are 
symmetric. This symmetry is not essential for the operation of the electronic 
neuron. For the given capacitances, the effective input capacitance, defined in (19 1, is 
C; n = 0.599 Co, which is why we have chosen C out = 0.6 Co- 

We measure time in units of the inverse rate, r^ 1 . To have a consistent unit 
system, we require RqCq — T^ 1 . If we choose a time scale of T^ 1 = 10 _10 s, which 
leads to feasible tunneling times [27l [28], and a unit capacitance of 10 _19 F, which 
appears to be realistic [231 GUI 1311 132] , we obtain R = T^ 1 /C = 1 GO. This resistance 
is much larger than the quantum resistance Rk — h/e 2 — 2.58 x 10 4 O of an open 
channel. Moreover, the input and output resistances in our circuit are much larger 
than i?o, which may be problematic. However, larger capacitances or faster tunneling 
would allow us to use smaller resistances. Our unit of energy is the energy quantum of 
the vibrational mode, eo = hui v . eo = 30 meV is a reasonable order of magnitude |27) . 
Since we will see that typical voltage signals are on the order of a few times eo / e, this 
leads to voltages on the order of 0.1 V. As noted above, in the simulations we assume 
a constant input current I- m , given in units of er^ 1 , and a vanishing output current 
lout = 0. The largest input current we use is |/i n | = 6x 10 _4 er o " 1 , which small in 
natural units. We first equilibrate the system for a time i; n ;t = 5x 10 6 T^ 1 , and then 
record observables for t sim — 10 7 Tq 1 , unless noted otherwise. 

Results for the ON state are given in figure [3] Here, the input current is 
-fin = —6 x 10 -4 erjj" 1 . Figure ga) shows the change in the number of electrons 
in the source electrode of the FC molecule, ATpcs- Clearly, the electrons flow in 
avalanches, typical for the FC regime. Comparison with figure 3 of [8] shows that the 
auxiliary molecular transistors do not appreciably change this behavior. In particular, 
feedback from the current through the second auxiliary molecule does not choke off 
the avalanche- like transport. 

In figuregb) we plot the number N^ot of electrons in the quantum dot, relative to 
the neutral state. Negative electron numbers simply mean that the dot is positively 
charged. The rapid fluctuations are due to electrons tunneling through the barrier 
between the dot and ground. The electron number preferentially assumes the values 
—4 and —5 since these values tune the onsite potential Vdot close to zero, which is the 
value it would assume in equilibrium with only the ground contact. The avalanches 
transmitted through the FC molecule are clearly visible. Note that the dot charge 
stays relatively small so that a large organic molecule or small metallic cluster should 
be able to accomodate it. 

Figure |3jc) shows the input and output voltages, V\ n and V on t, respectively. V\ n 
contains a constant term resulting from the constant input current. Sharp spikes are 
evident in both voltages, showing that the active neuron indeed fires voltage spikes 
into both the input and the output line. 

The spiky voltage V ln is the gate voltage seen by the FC molecular transistor. 
As we will see below, Vi n reaches values that would, if applied continuously, switch 
off the FC molecule. Nevertheless, the transport through the FC molecule is hardly 
affected by these spikes, as figure [3] shows. To understand this, we focus on a single 
series of spikes, taken from the same data as in figure § Figures ga) and gb) 
show Npc,s an d Vi n , respectively, for a shorter time interval. In figure gc) we plot 
the occupation number n c i and the harmonic-oscillator quantum number q of the FC 
molecule. Evidently, q > 1 during most of the episode. Thus the FC molecular 
transistor stores energy in the vibrational mode during the avalanches, which allows 
it to stay active even though the gate voltage is strongly reduced in magnitude. We 
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t (r ] ) 

Figure 3. Simulation of the molecular neuron in the ON state for the input 
current 7i n = —6 X 10 — ^eTg. The other parameters are given in the text, (a) 
Change in the number of electrons in the source electrode of the FC molecule, 
JVfc, s - (b) Number N^ ot °f electrons in the quantum dot. (c) Voltages V[ n and 
V ou t measured at the input and output contacts of the neuron, respectively. 

see that the memristive properties of the FC molecule are important: The current 
flowing through it does not only depend on the instantaneous gate and bias voltages 
but also on its history, here realized by the vibrational quantum number q |33j . 

Results for the OFF state are shown in figures |5|>,)-(c) for 7 in = 0. The same 
quantities as for the ON state in figure [3] are plotted. The FC molecule is not 
transmitting. The quantum dot only shows the fluctuating occupation number due 
to tunneling through the barrier, which is not sufficient to switch on the auxiliary 
molecular transistors. We also require the molecular neuron to show threshold 
behavior — it should not fire spikes for input currents below a certain threshold. Figures 
[5jd)-(f) show that the device can indeed take a finite input current without becoming 
active. Here, we have used L ln — —3 eT^ 1 , half the value for the active state in figurepl 

For neural functionality, it is desirable that the switching between OFF and ON 
states takes place over a narrow range of input currents. To study this onset, we 
plot in figure [6] the total number iV ou t of electrons inserted into the drain electrode 
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Figure 4. Details of the simulation presented in figure [3] for a shorter time 
interval, (a) Change in the number of electrons in the source electrode of the 
FC molecule, TVpc s- (b) Voltage V\ n measured at the input contact of the 
neuron. This voltage acts as the gate voltage of the FC molecular transistor, 
(c) Occupation number n,j and harmonic-oscillator quantum number q of the FC 
molecule. 



of the second auxiliary molecule, as a function of the input current A longer 
simulation time of t B i m = 5 x lO 7 ! 1 ^ 1 has been used. There is considerable noise, 
in particular for intermediate currents, since the FC molecule transmits electrons in 
avalanches. Between avalanches, the FC molecule returns to the vibrational ground 
state, q = 0, and thus does not retain any memory of the previous avalanche. Figure 
[6] essentially shows the shot noise of independent avalanches. Analysis of the time 
series (not shown) indicates that the main effect of tuning 7j n is to change the delay 
times between avalanches, not so much their duration. In any case, it is clear that the 
neuron is inactive over a wide range of input currents. 
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Figure 5. (a)-(c) Simulation of the molecular neuron in the OFF state for 
vanishing input current, 7; n = 0. The other parameters are the same as for 
figure [3] (a) Change in the number of electrons in the source electrode of 
the FC molecule, NpQ g . (b) Number N^ ot of electrons in the quantum dot. 
(c) Voltages V5 n and V ut measured at the input and output contacts of the 
neuron, respectively, (d)-(f) The same quantities obtained for the input current 
I ln = —3 X 10 — 4 er^" 1 . The other parameters are the same as for figure [3] The 
neuron is still inactive. 
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Figure 6. Number iVout of electrons inserted into the drain electrode of the 
second auxiliary molecule during a simulation time of t s i m = 5 X 10 7 r^" 1 , as a 
function of the input current Ii n . 
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5. Summary and conclusions 

We have proposed a design for a molecular-electronics realization of a neuron. The 
critical component of this artificial neuron is a molecular transistor with strong 
electron- vibration coupling, which is tuned to the FC-blockade regime when the neuron 
is active. In this regime, electrons flow in avalanches. The charge dumped by an 
avalanche into a small quantum dot or large molecule is used to activate two auxiliary 
molecular transistors, which lead to voltage pulses traveling along the output line 
and, importantly, also back up the input line. Employing Monte Carlo simulations for 
the dynamics of the circuit within the sequential-tunneling approximation, we have 
demonstrated that the device shows the desired neural behavior. Such a system can be 
fabricated with experimental capabilities available now or in the near future. It can be 
used for purposes other than the one discussed here. For instance, with an appropriate 
choice of system parameters it can operate as an amplifier. Most importantly, it can be 
used for dynamic storage and transfer of information in complex and highly integrated 
artificial networks that mimic the behavior of biological neural systems. 
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